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CN . Abstract 

I We present complete solutions of singlet and non-singlet Altarelli-Parisi (AP) evo- 

lution equations in leading order (LO) at low-x. We obtain t-evolutions of proton and 
0^ , neutron structure functions and x-evolutions of deuteron structure functions at low-x 

from AP evolution equations. The results of t-evolutions are compared with HERA 
low-x and low-Q^ data and those of x-evolutions are compared with NMC low-x and 
. low-Q^ data. Also we compare our results with those of general solutions of AP evolu- 

^ I tion equations. 
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1. Introduction: 

1^ • The Altarelli-Parisi (AP) evolution equations [1-4] are fundamental tools to study the 

t(=ln((5^/A^)) and x-evolutions of structure functions, where x and are Bjorken scaling 
^ ! and four momenta transfer in a deep inelastic scattering (DIS) process [5] respectively and 
^1 A is the QCD cut-off parameter. On the otherhand, the study of structure functions at 
low-x has become topical in view [6] of high energy collider and supercollider experiments 
[7]. Solutions of AP evolution equations give quark and gluon structure functions which 
produce ultimately proton, neutron and deuteron structure functions. Though numerical 
solutions are available in the literature [8], the explorations of the possibility of obtaining 
analytical solutions of AP evoluation equations are always interesting. In this connection, 
general solutions of AP evolution equations at low-x in leading order have already been ob- 
tained by applying Taylor expansion method [9] and t-evolutions [10] and x- evolutions [11] 
of structure functions have been presented. 

The present paper reports complete solutions of AP evolution equations in leading order 
at low-x and calculation of t and x-evolutions for singlet and non-singlet structure func- 
tions and hence t- evolutions of proton and neutron structure functions and x-evolutions of 
deuteron structure functions. In calculating structure functions, input data points have been 
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taken from experimental data directly unlike the usual practice of using an input distribu- 
tion function introduced by hand. Results of proton and neutron structure functions are 
compared with the HERA low-x low-Q^ data and those of deuteron structure functions are 
compared with the NMC low-x low-Q^ data. Comparisons are also made with the results of 
general solutions [10-12] of AP evolution equations already obtained. In the present paper, 
section 1 is the introduction. In section 2 necessary theory has been discussed. Section 3 
gives results and discussions. 



2. Theory: 

Though the basic theory has been discussed elsewhere [10-12], the essential steps of the 
theory has been presented here for clarity. The AP evolution equations for singlet and 
non-singlet structure functions in the standard forms are [13] 

^^^^^ - ^[{3 + 4/n(l - x)}Fi{x, t) + /f (x, t) + (x, t)] = (1) 

and 



- ^[{3 + Un{l - x)}F^^{x, t) + 7^^(x, t)\ = 0, (2) 



dt t 



where. 



js ^ 2 r ^{(1 + w^)Fi{^, t) - 2Fi{x, t)}, (3) 

Jx 1 — W W 



II = [\w' + (1 - wr}Gi-, t)dw (4) 



and 



/^^ = 2/ -^{{l + w^)F,^'{^,t)-2F,^'ix,t)}. (5) 

Here, t = ln{Q'^/A'^) and Af ^ 4/(33 - 2A^/), Nf being the number of flavours and A is the 
QCD cut off parameter. 

Using Taylor expansion method [9] and neglecting higher order terms as discussed in 
our earlier works [10-12], G{^,t), F2{^,t) and F^^{^,t) can be approximated for low-x as 

g(£*).g(.,«)+.|:.'^, (6) 



Fi(-,t) c^Fi{x,t)+xf2u 



^,dFi{x,t) 



dx 



(7) 



2 



and 



where x=l-w. 



Using equations (6), (7) and (8) in equations (3), (4) and (5) and performing u- integrations 
we get 

/f = -[(1 -x){x + ^)]Fl{x, t) + [2xln{-) + x{l - x^)]^^, (9) 



ll = ^/[^(l -x)i2-x + 2x'')G(x, t) + {~x(l -x)(5-Ax + 2x') + ^xZn(-)}^^^](10) 
2 / V / 2 ' ' ' 2 X ox 

and 

jNS ^ _ ^)(^ + 3)]F2^^(a;, t) + [2a;/n(i) + x{l - x^ /^^f'^^ (11) 
Now using equations (9) and (10) in equation (1) we have, 

- ^IA(.)FH., t) + B(.f-Mm. + CWG(.. *) + C = 0. (12) 

Let us assume for simphcity 

G{x,t)^K{x)Fi{x,t), (13) 

where K(x) is a function of x. Our earlier assumption as more crude, K=constant. This 
assumption gives t-evolutions of proton structure functions [8, 10, 12] which agree well with 
data. But it can not produce x-evolutions of deuteron structure functions [11] which agree 
well with data, especially in high-Q^ low-x region. Now equation (12) gives 

^?ffiM-^[LWf/.AfW^]^0. (14) 



where. 



A{x) = 3 + Aln{l - x) - (1 - x)(3 + x), 
Bix) =x(l-x2) + 2xZn(-), 

X 

C{x) = Nf{l -x){2-x + 2x^), 
D{x) = -\Nf{l - x)(5 - 4x + 2^2) + hn{-) 

L{x) = A{x) + K{x)C{x) + Dixf^^""^ 



dx 



and 

M{x) = B{x) + K{x)D{x). 
Secondly using equation (11) in equation (2) we have 

M;M_^[P(,.,,^«.^g(,.,)^?£rM, = 0. (15) 

where, 

P{x) = 3 + 4Zn(l -x)-{l-x){x + 3) 

and 

Q{x) = x{l — x"^) — 2xlnx. 
The general solutions of equation (14) is 

F(u, y) = 

where F is an arbitrary function and 

Uix,t,Fi)^Ci 

and 

V{x,t,Fi)^C2 



form a solution of equation 

dx dt dF^{x,t) 



(16) 



AfM{x) -t -AfL{x)Fi{x,ty 
Solving equation (16) we obtain, 

U(.,t,Fi)=te.plj-Jj^^d.] 

and 

V{x,t,Fi) = Fi{x,t) expij ^^dx]. 

If U and V are two independent solutions of equation (16) and if a and/3 are arbitrary 
constants, then 

V^aU + p 
is called a complete solution of equation (16). 



The complete solution 



Fiix, t) expiJ ^dx] = at exp[j- J jj^dx] + P 



M{x) ' "^'AfJ M{x) 
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is a two-parameter family of planes. The one parameter family determined by taking (3 — 
has equation 

/Tj( X) 1 /" X 

j^dx] = at exp[— J j^dx] + a\ (17) 

Differentiating equation (17) with respect to a, we get 

1 r 1 /■ 1 , 1 

a — —-t exp\—— I -———dx\. 

2 ^^Af J M{x) ^ 

Putting the value of a in equation (17), we obtain 

F/(..0 = -i«^..ri/(^-||)dx]. (18) 

Now, defining 

a,t t — to, where to — ln{Ql/ h?) at any lower value Q — Qo, we get from equation (18) 

Fi{x,t)^Fi{x,to){U^ (19) 
Jo 

which gives the t-evoluation of singlet structure function F2{x,t). 
Proceeding exactly in the same way, and defining 

we get for non-singlet structure function 

F,^^(x,^) = F,^^(x,^o)(f)^ (20) 

which gives the t-evolution of non-singlet structure function F^^{x,t). 
Again defining. 



\x=xo : 



we obtain from equation (18) 

Fi(., t) = t) eMjjjj^^ - (21) 



which gives the x-evoluation of singlet structure function F2{x,t). Similarly defining, 

we get 

Fi!'{x,t) = F^'{xo^) expiH-^^ - (22) 

Jxo AfQ{x) Q{x) 

which determines the x-evoluation of non-singlet structure function F^^{x,t). Deutron, 
proton and neutron structure functions measured in deep inelastic electro-production can be 
written in terms of singlet and non-singlet quark distribution functions in leading order as 

Fi{x,t)^^Fi{x,t), (23) 



Fiix, t) = ^Fiix, t) + ^Fi^H^, t) (24) 

and 

F-{x,t) = ^Fi{x,t) - ^Fi''{x,t). (25) 

Now using equations (19) and (21) in equation (23) we will get t and x-evolution of deuteron 
structure function ^^(x, i) at low-x as 

F,\x,t) ^ F^{x,to){lr (26) 



and 

Fi(., t) = t) -Pl£(^ - ^^Wl (27) 

where, the input functions are 



F^{xM)^\Fi{xM) 

and 

F^{xo,t) = ^-F^{x^,t). 

The corresponding results for general solutions of AP evolution equations obtained earlier 
[10-12] are 

F^{x,t) = F^{xM){h (28) 
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and 



fH., t) = F^(.„, t) eMjjjJ^^ - ^^)d4 (29) 

These have been obtained by taking arbitrary hnear combinations of U and V of arbitrary 
function F(U,V) — as solutions of equqations (14) and (15) for singlet and non-singlet 
structure functions respectively. 



Similarly using equations (19) and (20) in equations (24) and 25) we get the t-evolutions 
of proton and neutron structure functions at low-x as 

Fi{x,t) = Fi{x,to){lr (30) 



and 



F^{x,t) = F^{x,to){h\ (31) 

where, the input functions are 

Fiix,to) = ^Fi{x,to) + ^F2^^(x,to) 

and 

Fn^,to) = ^Fi{x,to) ^ ^Ff^(x,to). 
The corresponding results for general solutions of AP evolution equations are 

FUx,t) ^ FUx,to){l) (32) 

and 

F^{x,t) ^ F^{x,to){U. (33) 
to 

But the x-evoluations of proton and neutron structure functions like those of deuteron struc- 
ture function is not possible by this methodology, because to extract the x-evoluation of 

proton and neutron structure functions wc arc to put equations (21) and (22) in equations 
(24) and (25). But as the functions inside the integral sign of equations (21) and (22) are 
different, we need to separate the input functions F2{xo,t) and F^'^{xQ,t) from the data 
points to extract the x-evolutions of the proton and neutron structure functions, which is 
not possible. 



3. Results and Discussion: 

In the present paper, we compare our results of t-evolutions of proton and neutron struc- 
ture functions from equations (30) and (31) respectively with the HERA low-x, low-Q^ data 
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[14]. Here proton structure functions (x, Q^, z) measured in the range 2 < < 50 GeV'^, 
0.73 < z < 0.88 and neutron structure functions F^ix^Q"^ . z) measured in the range 
2 < g2 < 50 GeV^, 0.3 < 2 < 0.9 have been used. Moreover here Pt < 200 Mev, 
where Pt is the transverse momentum of the final state baryon and z = 1 — q-{p — p')/ (q-p), 
where p, q are the four momenta of the incident proton and the exchanged vector boson 
couphng to the positron and p' is the four-momentum of the final state baryon. 

In figures l(a-d) we present our results of t-evolutions of proton structure functions Ff 
(solid hnes) for the representative values of x given in the figures. Data points at lowest-Q^ 
values in the figures are taken as input to test the evolution equation (30). Agreement is 
found to be excellent. In the same figures we also plot the results of t-evolutions of proton 
structure functions (dashed lines) for the general solutions from equation (32) of AP 
evolution equations. We observe that results of complete solutions are better than those of 
general solutions. 

In figures 2(a-d) we present our results of t-evolutions of neutron structure functions F^ 
(solid lines) for the representative values of x given in the figures. Data points at lowest-Q^ 
values in the figures are taken as input to test the evolution equation (31). Agreement is 
found to be excellent. In the same figures we also plot the results of t-evolutions of neutron 
structure functions F|^ (dashed lines) for the general solutions from equation (33) of AP evo- 
lution equations. We observe that in this case also results of complete solutions are better 
than those of general solutions. 

For a quantitative analysis of x-distributions of structure functions, we evaluate the 
integrals that occured in equation (27) and present the results in figures 3(a-d) for Nf = 4 
and representative values of given in each figure, and compare them with NMC deuteron 
low-x low-Q^ data [15]. In each figure the data point for x- value just below 0.1 has been taken 
as input F2 (xo, t). If we take K(x)=constant, agreement of the results with experimental data 
becomes poor as before [11] especially in high-Q^ lower-x region. We, therefore, consider here 
K(x)=a.exp(-bx) where a=l and b=l for simplicity which gives better results. But explicit 
form of K(x) can actually be obtained only by solving coupled AP evolution equations for 
singlet and gluon structure functions, and works are going on in this regard. 

Traditionally the AP equations provide a means of calculating the manner in which the 
parton distributions change at fixed x as varies. This change comes about because of the 
various types of parton branching emission processes and the x-distributions are modified as 
the initial momentum is shared among the various daughter partons. However the exact rate 
of modifications of x-distributions at fixed Q'^ can not be obtained from the AP equations 
since it depends not only on the initial x but also on the rate of change of parton distribu- 
tions with respect to x, dP'F[x)/dx^ {u—1 to 00), up to infinite order. Physically this implies 
that at high-x, the parton has a large momentum fraction at its disposal and as a result it 
radiates partons including gluons in innumerable ways, some of them involving complicated 
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QCD mechanisms. However for low-x, many of the radiation processes will cease to occur 
due to momentum constraints and the x-evoluations get simplified. It is then possible to 
visualise a situation in which the modification of the x-distribution simply depends on its 
initial value and its first derivative. In this simplified situation, the AP equations give infor- 
mation on the shapes of the x-distribution as demonstrated in this paper. The clearer testing 
of our results of x-evolution is actually the equation (22) which is free from the additional 
assumption equation (13). But non-singlet data is not sufficiently available in low-x to test 
our result. It is observed as a whole that the results of complete solutions of AP evolution 
equations are improved than those of general solutions, especially in t-evolution calculations. 
Of course, this is a leading order calculation. Its natural improvement will be the calculation 
considering next-to-leading order terms and we plan to do so in our subsequent works. 
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Figure 1: In figures l(a-d) t-evolutions of proton structure functions F2 (solid lines) for the 
representative values of x given in the figures. Data points at lowest-Q^ values in the figures 
are taken as input to test the evolution equation (30). To test the evolution equation (30). 
In the same figures we also plot the results of t-evolutions of proton structure functions 
(dashed hues) for the general solutions from equation (32) of AP evolution equations. 
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Figure 2: In figures 2(a-d) t-evolutions of neutron structure functions F2 (solid lines) for 
the representative values of x given in the figures. Data points at lowest-Q^ values in the 
figures are taken as input to test the evolution equation (31). In the same figures we also plot 
the results of t-evolutions of neutron structure functions F2 (dashed fines) for the general 
solutions from equation (33) of AP evolution equations. 
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Fig'3(a), Q^=0.75 GeV^ 
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Figure 3: In figures 3(a-d) x-distributions of deuteron structure functions (solid lins) for 
Nf — 4 and the representative values of given in each figure and compared them with 
NMC deuteron low-x low Q"^ data. In each figure the data point for x-value just below 0.1 
has been taken as input F2{xo,t). 
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